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ABSTRACT 

The  problem  of  association  and  fusion  of  radar  and  infra-red  search  and  track 
(IRST)  sensor  reports  is  not  straightforward,  especially  because  IRST  pro¬ 
vides  only  angular  measurements  while  radar  provides  range  and  range-rate 
measurements  in  addition  to  azimuth  and  elevation  measurements.  In  this  re¬ 
port,  simulation  results  show  that  a  centralised  extended  Kalman  filter  tracker 
is  a  solution  that  can  capitalise  on  the  higher  angular  accuracy  of  the  IRST 
sensor  to  provide  improved  track  accuracy  performance. 
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Tracking  Anti-Ship  Missiles  Using  Radar  and  Infra-Red 
Search  and  Track:  Track  Error  Performance 


EXECUTIVE  SUMMARY 

Modern  anti-ship  missiles  represent  a  serious  threat  to  maritime  assets.  They  are 
designed  to  fly  close  to  the  sea  surface  and  can  perform  high-acceleration  terminal  ma¬ 
noeuvres,  making  them  difficult  to  detect  and  track  in  time  for  effective  engagement  by 
shipboard  weapons  systems. 

Following  two  decades  of  research  and  development  in  ship  self  defence  systems,  one  of 
the  promising  aspects  that  has  emerged  is  the  combination  of  radar  and  infra-red  search 
and  track  (IRST)  sensors  to  improve  the  detection  and  tracking  of  anti-ship  missiles. 

This  is  recognised  in  Project  SEA  1448,  which  is  required  to  deliver  enhanced  ship 
self  defence  against  modern  anti-ship  missiles  for  the  Royal  Australian  Navy’s  ANZAC 
frigates.  As  part  of  this  project,  IRST  sensors  will  be  integrated  into  an  upgraded  combat 
system. 

In  practice,  implementation  of  effective  operational  data  fusion  systems  is  far  from 
simple  and  fusion  of  sensor  data  may  actually  produce  worse  results. 

In  this  report,  Monte  Carlo  simulations  are  employed  to  quantify  the  track  error  achiev¬ 
able  using  standard  tracking  algorithms  (namely,  global  nearest  neighbour  data  association 
and  an  extended  Kalman  filter)  with  and  without  IRST  when  tracking  non-manoeuvring 
and  manoeuvring  missile  targets.  The  aim  is  to  determine  whether  combining  radar  with 
IRST  provides  significant  benefits  in  reducing  track  error. 

Three  missile  trajectories  are  considered: 

Straight  The  missile  flies  in  a  straight  line  at  constant  speed  and  height  towards  the  ship. 

Weave  The  missile  flies  in  a  sinusoidal  line  at  constant  height  towards  the  ship. 

Dive  The  missile  flies  in  a  straight  line  towards  the  ship  and  then  rapidly  climbs  to  a  new 
height  in  preparation  for  a  dive  attack. 

Track  error  performance  is  measured  in  terms  of  90%  azimuth,  elevation  and  range  error 
levels  from  100  Monte  Carlo  simulation  runs  for  each  of  the  three  missile  trajectories. 

The  simulation  results  confirm  that  the  incorporation  of  IRST  can  lead  to  improved 
track  accuracy  performance.  The  track  azimuth  and  elevation  errors  drop  to  the  IRST 
measurement  error  level,  while  the  track  range  error  does  not  change  significantly.  This  is 
to  be  expected  since  IRST  measures  azimuth  and  elevation  but  not  range. 
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1  Introduction 

Modern  anti-ship  missiles  represent  a  serious  threat  to  maritime  assets.  They  are 
designed  to  fly  close  to  the  sea  surface  and  can  perform  high- acceleration  terminal  ma¬ 
noeuvres,  making  them  difficult  to  detect  and  track  in  time  for  effective  engagement  by 
shipboard  weapons  systems. 

Following  two  decades  of  research  and  development  in  ship  self  defence  systems,  one  of 
the  promising  aspects  that  has  emerged  is  the  combination  of  radar  and  infra-red  search 
and  track  (IRST)  sensors  to  improve  the  detection  and  tracking  of  anti-ship  missiles. 

This  is  recognised  in  Project  SEA  1448,  which  is  required  to  deliver  enhanced  ship 
self  defence  against  modern  anti-ship  missiles  for  the  Royal  Australian  Navy’s  ANZAC 
frigates.  As  part  of  this  project,  IRST  sensors  will  be  integrated  into  an  upgraded  combat 
system. 

Some  of  the  benefits  of  combining  radar  with  IRST  are  [Missirian  &  Ducruet  1997, 
Misanin  1995] 

Increased  detection  performance  Radar  performance  is  degraded  at  low  elevations 
due  to  surface  clutter,  multipath  propagation  and  unfavourable  evaporative  ducting 
conditions.  IRST  operates  at  different  frequency  bands  and  provides  additional 
means  of  target  detection  via  aerodynamic  skin  heating  and  engine  exhaust  plumes. 

Reduced  threat  alert  time  For  an  incoming  missile,  skimming  along  the  surface  of  the 
sea,  the  time  between  when  the  missile  becomes  visible  on  the  radar  horizon  and 
when  it  reaches  its  target  is  very  short.  The  time  can  be  less  than  30  seconds  for  some 
types  of  anti-ship  missile  [Horman,  Stapleton,  Hepfer,  Headley  &  Stapleton  1996]. 
An  IRST,  performing  dedicated  horizon  search,  can  provide  earlier  target  detection 
and  cue  a  radar  to  the  target  for  more  observations. 

Reduced  vulnerability  to  ECM  Radar  emissions  are  restricted  when  covertness  is  re¬ 
quired,  or  to  minimise  susceptibility  to  jamming.  In  such  circumstances,  an  IRST 
may  become  the  primary  sensor  to  guard  against  anti-ship  missiles. 

Another  potential  benefit  of  combining  radar  and  IRST  sensors  is  improved  track 
accuracy  which  is  important  for  sensor  cueing  and  weapons  control  functions.  By  appro¬ 
priate  means  of  multisensor  data  fusion,  the  radar’s  accurate  range  measurements  (but 
not-so-accurate  angle  measurements)  can  be  combined  with  the  IRST’s  accurate  angle 
measurements  to  provide  a  fused  estimate  of  the  target’s  location  that  is  of  reduced  error 
[Hall  &  Llinas  1997]. 

In  practice,  implementation  of  effective  operational  data  fusion  systems  is  far  from 
simple  and  fusion  of  sensor  data  may  actually  produce  worse  results  [Hall  &  Llinas  1997]. 

Other  authors  have  reported  track  error  performance  using  IRST  alone  [Maltese  2001], 
enhanced  angular  track  accuracy  achievable  by  fusing  IRST  and  ESM  sensors  [Maltese  & 
Lucas  2000],  the  effect  of  the  choice  of  coordinate  system  on  track  error  performance  of 
radar  and  IRST  fusion  [Simard  &  Begin  1993],  and  significant  reduction  in  track  error 
by  applying  Multiple  Hypothesis  Tracking/Interacting  Multiple  Models  (MHT/IMM)  on 
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radar  and  IRST  data  [Blackman,  Dempster  &  Roszkowski  1997,  Dempster,  Blackman, 
Roszkowski  &  Sasaki  1998]. 

In  this  report,  Monte  Carlo  simulations  are  employed  to  quantify  the  track  error 
achievable  using  standard  tracking  algorithms  with  and  without  IRST  when  tracking 
non-manoeuvring  and  manoeuvring  missile  targets.  The  aim  is  to  determine  whether 
combining  radar  with  IRST  provides  significant  benefits  in  reducing  track  error. 

Only  generic  radar  and  IRST  sensor  models  are  used  and  no  attempt  has  been  made 
to  model  environmental  conditions.  Simulated  anti-ship  missile  trajectories  include  weave 
and  dive.  Further  details  of  the  simulated  missile  trajectories,  sensor  models  and  tracking 
algorithms  are  given  in  Sections  2,  3  and  4  respectively.  The  error  performance  statistic  is 
defined  in  Section  5.  The  simulation  results  are  given  and  discussed  in  Section  6.  Section  7 
provides  the  conclusions. 


2  Missile  Trajectories 

Three  possible  missile  trajectories  are  considered: 

Straight  The  missile  flies  in  a  straight  line  at  constant  speed  and  height  towards  the  ship. 

Weave  The  missile  flies  in  a  sinusoidal  line  at  constant  height  towards  the  ship. 

Dive  The  missile  flies  in  a  straight  line  towards  the  ship  and  then  rapidly  climbs  to  a  new 
height  in  preparation  for  a  dive  attack. 

The  ship  is  located  at  the  origin.  Each  of  the  simulated  trajectories  lasts  for  27  seconds. 
Specific  details  are  given  in  the  subsections  that  now  follow. 

2.1  Straight 

Specific  parameters,  relative  to  the  ship  origin,  are: 

Bearing  90° 

Start  Range  19  km 
Speed  700  m/s 
Height  10  m 

2.2  Weave 

Specific  parameters,  relative  to  the  ship  origin,  are: 

Bearing  0° 
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Figure  1:  Weave  missile  speed  and  acceleration 


Start  Range  19  km 
Height  10  m 
Weave  Period  5  s 
Weave  Amplitude  89  m 

Figure  1  gives  the  instantaneous  speed  and  acceleration  of  the  weave  missile.  It  can  be 
seen  that  acceleration  reaches  just  above  15  g  (where  1  g  equals  9.81  m/s2). 

2.3  Dive 

The  dive  trajectory  is  based  on  the  bunt  manoeuvre  scenario  of  Maltese  (2001).  Fig¬ 
ure  2  gives  the  trajectory  simulated  for  27  seconds  only.  Most  of  the  dive  portion  occurs 
after  27  seconds  and  is  therefore  not  apparent  in  the  figure.  Figure  3  gives  the  instan¬ 
taneous  speed  and  acceleration  of  the  dive  missile.  It  can  be  seen  that  acceleration  can 
exceed  20  g  during  the  height  increase  and  30  g  as  the  missile  dives. 


3  Sensor  Models 

The  radar  and  IRST  sensor  are  characterised  by  range-dependent  “probability  of  target 
detection”  curves  for  a  constant  probability  of  false  alarm.  Measurement  accuracy  is 
modelled  by  the  Gaussian  distribution  for  the  sensor  noise.  Both  sensors  rotate  360°  and 
have  the  scan  times  as  given  in  Table  1.  Note  that  two  surveillance  configurations  are 
modelled,  radar-only  and  radar-plus-IRST  configurations. 
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Figure  2:  Dive  trajectory 


Figure  3:  Dive  missile  speed  and  acceleration 


Table  1:  Sensor  scan  times 


Configuration 

Radar 

IRST 

Radar-only 

Radar-plus-IRST 

5  s 

5  s 

1  s 
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Table  2:  Radar  measurement  accuracy 


Measurement 

Error  standard  deviation 

Range 

60  m 

Azimuth 

0.5° 

Elevation 

1° 

Range-rate 

0.7  m/s 

Table  3:  Radar  field  of  view 


Azimuth 

360° 

Elevation 

O 

o 

co 

1 

o 

O 

Range 

1-450  km 

Range-rate 

0-1000  m/s 

The  sensors  are  ideal  in  that  they  have  negligible  registration  error  and  they  resolve 
targets  fully. 

Note  that  environmental  effects  such  as  clutter,  multipath  propagation  and  evaporative 
ducts  are  not  modelled.  These  effects  could  degrade  sensor  detection  performance  in 
practice  and  reduce  tracking  performance  as  a  result.  Future  work  could  look  at  the 
impact  of  these  effects  on  tracking  performance. 


3.1  Radar 

For  rapidly-approaching  targets,  it  is  assumed  that  the  signal-to-noise  ratio  (SNR)  for 
targets  is  not  attenuated  by  the  radar  processor’s  Doppler  notch. 

For  a  10-5  probability  of  false  alarm,  the  probability  of  detection  is  1.0  for  a  Swerling  1 
target  at  ranges  less  than  19  km,  according  to  textbook  formulae  [Blackman  &  Popoli  1999, 
sec.  2.2]. 

If  the  radar  detects  a  target,  it  returns  noisy  measurements  of  the  azimuth,  elevation, 
range  and  range-rate  of  the  target.  Table  2  gives  the  standard  deviations  chosen  for  the 
measurement  noise.  False  alarm  measurements  are  uniformly  distributed  over  the  radar 
field  of  view.  The  radar  field  of  view  is  given  in  Table  3. 


3.2  IRST 

The  IRST  sensor  is  an  array  of  640  by  512  detector  elements  which  captures  the  infra¬ 
red  scene  once  every  1.5  ms.  For  a  target  at  range  i?,  the  IRST  measures  a  signal  with 
amplitude  that  is  proportional  to  R~ 2  [Blackman  &  Popoli  1999,  sec.  2.3.7].  This  signal 
is  corrupted  by  Gaussian  noise.  A  threshold  is  set  to  maintain  a  false  alarm  rate  of  not 
more  than  1.0  per  hour  over  the  entire  array  of  detector  elements.  Given  this  threshold, 
the  signal  strength  from  the  target  is  such  that,  at  R  —  11  km,  the  probability  of  detection 
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Figure  4:  The  probability  of  target  detection  in  every  scan  of  the  IRST 


Table  4:  IRST  measurement  accuracy 


Measurement 

Error  standard  deviation 

Azimuth 

Elevation 

0.5  mrad 

0.5  mrad 

is  0.9.  Figure  4  depicts  the  resulting  single-scan  target  detection  probability  curve  for  the 
IRST. 

If  a  target  is  detected,  the  IRST  measures  the  azimuth  and  elevation  of  the  target. 
Table  4  gives  the  values  chosen  for  the  IRST  measurement  noise  standard  deviations.  For 
a  false  alarm,  the  measured  azimuth  and  elevation  are  uniformly  distributed  over  the  IRST 
field  of  view.  The  IRST  field  of  view  is  given  in  Table  5. 


Table  5:  IRST  field  of  view 


Azimuth 

360° 

Elevation 

0°-10° 

Range 

0.1-30  km 
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4  Target  Tracking 

Target  tracking  is  performed  in  rectangular  coordinates.  The  track  state  is  composed 
of  position  and  velocity  estimates  of  the  target. 


Reports  from  radar  and  IRST  sensors  are  sent  directly  to  a  central  tracker  that  is 
based  on  the  Extended  Kalman  Filter  (EKF). 

The  tracker  comprises  the  following  processes: 

•  Data  association 

•  Track  initiation 

•  Track  update 

•  Track  deletion 

Details  of  these  are  given  below. 

4.1  Data  Association 

Before  a  track  can  be  updated  by  a  report,  the  report  must  first  be  associated  with 
the  track.  This  process  generally  involves  calculating  a  correlation  statistic  for  each  track- 
report  pair  and  assigning  reports  to  tracks  based  on  these  correlation  statistics. 

The  gating  of  reports  also  forms  part  of  the  process.  Gating  is  applied  to  determine 
which  reports  are  reports  of  new  targets,  i.e.,  targets  which  have  not  been  detected  pre¬ 
viously.  For  reports  of  new  targets,  new  tracks  are  created  through  the  process  of  track 
initiation. 

A  basic  assignment  method  is  the  Global  Nearest  Neighbour  (GNN)  method  which 
considers  optimal  assignment  of  reports  to  tracks  in  a  fixed  frame  of  time.  A  more  so¬ 
phisticated  method  is  the  Multiple  Hypothesis  Tracking  (MHT)  method  which  considers 
optimal  assignment  over  multiple  frames  of  time.  MHT  is  computationally  intensive  but 
is  useful  for  tracking  closely-spaced  targets.  This  is  where  single-frame  methods  have 
difficulty  making  correct  assignments.  Further  details  of  these  methods  can  be  found  in 
standard  textbooks. 

In  this  report,  the  tracker  applies  GNN  logic  with  a  frame  time  of  0.5  seconds.  MHT 
was  not  considered  to  be  necessary  for  the  single  target  scenarios  of  interest. 
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The  correlation  statistic  for  a  track  with  state  x  and  covariance  P,  and  a  report  with 
measurement  z  and  measurement  covariance  R  is  given  by 

G  -  (z  -  h(x))'(HPH'  +  R)-1(z  -  h(x))  (2) 

where  h  is  the  measurement  function,  H  is  the  Jacobian  of  the  measurement  function, 
and  G  is  the  gating  threshold.  Note  that  the  track  state  and  covariance  are  quantities 
predicted  from  the  time  of  the  most  recent  track  update  to  the  time  of  the  report  being 
correlated.  More  details  are  given  in  Section  4.3  on  the  track  update  process. 

The  gating  threshold  G  is  an  adjustable  parameter  that  discriminates  reports  of  new 
targets  from  those  of  existing  targets.  It  is  set  to  105  in  the  EKF  tracker  for  gating  both 
radar  or  IRST  reports. 

The  Auction  algorithm  is  used  to  assign  the  reports  to  tracks.  This  algorithm  is 
described  in  Section  6.5.1  of  Blackman  Sz  Popoli. 


4.2  Track  Initiation 

Each  report  that  is  not  associated  with  a  track  after  the  data  association  process  is 
converted  from  spherical  to  rectangular  coordinates  to  provide  the  initial  track  position 
and  track  position  covariance.  To  minimise  the  biases  introduced  by  the  conversion,  the 
formulae  in  Appendix  A  are  applied. 

For  IRST  reports,  range  is  not  measured.  In  this  case,  the  conversion  assumes  a 
‘measured’  range  of  50  km  with  a  ‘measurement  error’  standard  deviation  of  50  km.  (As 
suggested  by  the  referee  of  this  report,  an  initial  range  of  30  km  might  be  a  more  appro¬ 
priate  assumption  given  that  the  IRST  field  of  view  for  range  is  0-30  km.) 

The  initial  track  velocity  is  set  to  zero  and  the  initial  track  velocity  standard  deviation 
is  set  to  267  m/s  in  each  coordinate. 

To  incorporate  the  range-rate  measurement  in  radar  reports,  a  single  EKF  update  is 
performed.  The  steps  are  given  in  Table  6. 

The  initial  track  state  and  covariance  represent  a  tentative  track  which  would  not 
be  output  to  an  operator  display  process  because  it  may  be  a  false  track.  Displaying 
such  tentative  tracks  would  lead  to  an  unacceptably  high  number  of  false  tracks.  This  is 
especially  true  in  heavy  clutter  environments  where  there  are  many  spurious  reports. 

The  track  becomes  an  established  track  after  it  is  updated  by  two  more  correlated 
reports.  This  track  is  then  ‘safe’  to  output  because  it  is  more  likely  to  represent  a  target 
of  interest  than  to  be  a  false  track. 


4.3  Track  Update 

When  an  existing  track  is  associated  with  a  report  by  the  data  association  process, 
it  is  updated  via  a  two-step  procedure.  The  first  step  of  this  procedure  is  to  predict  the 
track  to  the  time  of  the  report.  The  second  step  is  to  correct  the  predicted  track  using 
the  measured  data. 
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Table  6:  Steps  to  incorporate  a  range-rate  measurement  into  the  initial  track  state  and 
covariance 


Compute  Jacobian  matrix 

H  —  [(x  —  rx/r)/r  ( y  —  ry/r)/r  ( z  —  rz/r)/r  x/r 

y/r  z/r]  (3) 

where  (x,y,z)  and  (i,?/,i)  are  the  initial  track  state  position  and  velocity 

before  the  range-rate  update,  and 

r  —  \/  x1  +  y2  +  z2 

(4) 

r  —  {xx  +  yy  +  zz)/r 

(5) 

Compute  Kalman  gain 

K  =  PH'(HPH' +  er?)-1 

(6) 

where  P  is  the  initial  track  covariance  before  the  range-rate  update,  and 

df  is  the  range-rate  measurement  error  standard  deviation 
known. 

,  assumed  to  be 

Update  initial  track  state 

x  * —  x  K (z  —  r) 

(7) 

where  z  is  the  range-rate  measurement. 

Update  initial  track  covariance 

P  <-  P  -  KHP 

(8) 
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The  update  procedure  for  the  EKF  tracker  is  given  in  Table  7.  This  procedure  assumes 
a  white  noise  constant  velocity  target  [Blackman  &  Popoli  1999,  sec.  4.2.2].  For  radar 
reports,  it  converts  the  measured  range,  azimuth  and  elevation  to  rectangular  coordinates 
using  the  formulae  in  Appendix  A.  For  IRST  reports,  the  update  procedure  is  given  in 
Table  8. 

Note  that  the  process  noise  parameter  q  (in  (11)  of  Table  7)  is  set  empirically  to  obtain 
some  ability  for  the  EKF  tracker  to  track  through  target  manoeuvres  whilst  maintaining 
an  acceptable  track  error  performance  during  the  non-manoeuvring  sections  of  target 
trajectories. 

If  improved  performance  on  manoeuvring  targets  is  desired,  the  Interacting  Multiple 
Models  (IMM)  tracking  approach  is  recommended.  This  is  a  standard  approach  where  the 
target  motion  is  modelled  as  a  stochastic  process  that  can  switch  between  several  different 
dynamic  models.  For  example,  a  target  moving  at  constant  velocity  for  a  period  of  time 
and  then  switching  to  constant  acceleration  motion  can  be  modelled.  An  IMM  tracker  can 
provide  improved  track  accuracy  performance  at  increased  computational  cost  compared 
with  an  EKF  tracker. 


4.4  Track  Deletion 

Tracks  are  deleted  if  they  have  not  been  updated  for  more  than  10  seconds,  or  if  the 
predicted  track  error  variance  (from  the  diagonal  of  the  track  covariance  matrix  P)  exceeds 
105  m2  in  the  x  or  y  coordinate,  or  5000  m2  in  the  z  coordinate.  Under  these  conditions, 
tracks  are  said  to  be  lost  or  to  have  ‘gone  stale’. 


5  Performance  Measure 

In  this  report,  track  error  performance  is  measured  in  terms  of  90%  azimuth,  elevation 
and  range  error  levels  from  100  Monte  Carlo  simulation  runs  for  each  of  the  scenario 
described  in  Section  2. 

For  each  Monte  Carlo  simulation  run,  the  track  azimuth/elevation/range  error  (i.e., 
absolute  difference  between  the  track  and  true  values)  is  computed  and  sampled  at  every 
0.2  seconds.  If  there  are  multiple  tracks  at  a  given  sample  time,  the  smallest  track  error 
at  that  time  is  taken.  If  no  track  exists  at  a  given  sample  time,  no  error  sample  is  taken. 

Error  samples  collected  at  each  time  point  are  sorted  in  ascending  order.  The  90% 
error  level  bounds  90%  of  the  error  samples.  It  is  taken  as  the  Mth  sorted  error  sample, 
where  M  =  |_0.9(7V  +  1)J  and  N  denotes  the  number  of  samples. 


6  Results  and  Discussion 

Figures  5,  6  and  7  show  the  track  azimuth,  elevation  and  range  errors  for  straight, 
weave  and  dive  missile  trajectories  respectively.  Also  shown  are  the  radar  and  IRST 


10 


DSTO-TR-1863 


Table  7;  EKF  track  update  procedure 


Predict 


x  Fx 

P  <-  FPF'  +  Q 


(9) 

(10) 


where 


(ii) 


where  q  —  100,  and  T  is  the  time  increment  to  predict  the  track  to  the 
time  of  the  report.  I3  denotes  a  3-by-3  identity  matrix,  and  O3  denotes  a 
3-by-3  zero  matrix. 


Correct  For  radar  reports  which  have  range,  azimuth,  elevation  and  range-rate 
measurements: 


Compute  Jacobian  matrix 

1  0  0  0  0  0  " 

0  1  0  0  0  0 

0  0  1  0  0  0 

(. x  —  rx/r)/r  ( y  —  ry/r)/r  ( z  —  rz/r)/r  x/r  y/r  z/r 

. .  '  (12) 

where  ( x,y,z )  and  (rr,y,  i)  are  the  predicted  track  state  position  and 
velocity,  and 

r  =  \/  x2  +  y2  +  z2  (13) 

r  =  (xx  +  yy  +  zz)/r  (14) 


Compute  Kalman  gain 


K  =  PH(HPH'  +  R)1  (15) 

where  P  is  the  predicted  track  covariance,  and  R  is  the  measurement 
covariance. 

Update  predicted  track  state 


x^x  +  K(z  —  [xy  z  r]') 


(16) 


where  z  is  the  measurement  vector  comprising  of  the  3D  position  and 
measured  range-rate.  The  3D  position  is  in  rectangular  coordinates 
and  is  obtained  by  converting  the  spherical  measurements  (of  range, 
azimuth  and  elevation)  using  formulae  in  Appendix  A. 

Update  predicted  track  covariance 


P  <-  P  -  KHP 


(17) 
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Table  8:  EKF  track  update  procedure  for  IRST  reports 


Predict  the  track  as  in  Table  7  and  then  correct  the  predicted  track  as  follows: 

Compute  Jacobian  matrix 

-y/rg  x/rg  0  0  0  0 

xc  ye  rg/r 2  0  0  0 

where  (x,y,z)  is  the  predicted  track  state  position,  and 


H 


(18) 


r„  =  ^/x2  +  y2 


r  =  +  y2  +  z2 

c  =  -z/(r2rg) 


Compute  Kalman  gain 


K  =  PH'  HPH'  + 


Gn 


0  <r2 


-1 


(19) 

(20) 
(21) 


(22) 


where  P  is  the  predicted  track  covariance,  and  cra  and  ae  are  azimuth  and 
elevation  measurement  error  standard  deviations,  assumed  to  be  known. 


Update  predicted  track  state 


x  <-  x  +  K 


a  —  tan  1(y,x) 
e  —  tan_1(z,  rg) j 


(23) 


where  a  and  e  are  the  measured  azimuth  and  elevation.  Note  that  if 
| a  —  tan_1(?/,x)|  >  i r,  the  appropriate  multiple  of  27r  should  be  added  or 
subtracted  from  a  —  tan ~1(y,x)  to  obtain  a  value  between  —  tv  and  n  to 
replace  it  in  the  above  equation. 


Update  predicted  track  covariance 

P  <-  P  -  KHP 


(24) 
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- radar  measurement  error 

- IRST  measurement  error 

- radar-only  track  error 

— • —  radar+IRST  track  error 


Figure  5:  Track  error  for  the  straight  missile  trajectory 


measurement  errors  for  comparison.  These  errors  correspond  to  the  90%  quantile  of  the 
absolute  error  distribution  assuming  a  Gaussian  distribution  for  the  sensor  error  with  the 
standard  deviations  given  in  Table  2  and  Table  4. 

The  results  confirm  that  the  incorporation  of  IRST  can  lead  to  improved  track  accuracy 
performance.  The  track  azimuth  and  elevation  errors  drop  to  the  IRST  measurement  error 
level,  while  the  track  range  error  does  not  change  significantly.  This  is  to  be  expected  since 
IRST  measures  azimuth  and  elevation  but  not  range. 

Although  not  the  focus  of  this  report,  the  results  also  indicate  earlier  track  initiation 
due  to  the  faster  scan  rate  of  the  IRST  sensor.  This  is  apparent  in  the  figures.  The 
accuracy  performance  for  radar+IRST  tracks  is  available  after  at  most  7  seconds  into  the 
simulation  compared  with  10-11  seconds  for  radar-only  tracks. 


6.1  Straight 

The  results  in  Figure  5  demonstrate  the  value  of  sensor  fusion.  Track  angle  errors  are 
reduced  to  the  level  of  that  of  the  more  accurate  IRST  sensor,  while  track  range  error  is 
almost  unchanged. 


13 


DSTO-TR-1863 


- radar  measurement  error 

- IRST  measurement  error 

- radar-only  track  error 

— • —  radar+IRST  track  error 


Figure  6:  Track  error  for  the  weave  missile  trajectory 


6.2  Weave 

For  a  weaving  missile  that  is  manoeuvring  predominantly  in  the  azimuth  axis,  it  can 
be  seen  that  the  track  error  along  this  axis  is  elevated  compared  with  the  performance  for 
the  straight  trajectory.  There  is  some  slight  effect  of  the  weaving  on  radar+IRST  track 
range  error  which  is  noticeable  in  Figure  6. 


6.3  Dive 

In  the  case  of  the  dive  trajectory  where  the  missile  executes  rapid  changes  in  eleva¬ 
tion,  the  increased  scan  rate  of  the  IRST  sensor  helps  the  EKF  tracker  maintain  track 
on  the  target  and  keep  the  track  elevation  accuracy  close  to  the  accuracy  of  the  IRST 
measurement. 

Curiously,  the  results  for  the  dive  trajectory  do  not  exhibit  the  same  behaviour  com¬ 
pared  with  the  results  for  the  straight  and  weaving  trajectories  near  the  end  of  the  simu¬ 
lation.  That  is,  the  track  errors  for  the  dive  trajectory  do  not  shoot  off  the  graph  after  27 
seconds.  The  reason  for  this  is  that  the  missile  is  still  about  4  km  away  from  the  sensors  in 
this  scenario.  In  the  straight  and  weaving  missile  scenarios,  the  missile  is  only  a  hundred 
metres  or  so  away  after  27  seconds  have  elapsed.  At  such  a  close  range,  the  EKF  tracker 
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- radar  measurement  error 

- IRST  measurement  error 

- radar-only  track  error 

— • —  radar+IRST  track  error 


Figure  7;  Track  error  for  the  dive  missile  trajectory 
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cannot  cope  with  the  nonlinearity  of  the  measurements. 

In  practice,  targets  at  very  close  range  may  be  detected  in  more  than  one  sensor 
resolution  cell  and  sensor  reports  may  become  unreliable  or  unavailable.  It  is  hoped,  of 
course,  that  a  missile  never  gets  within  4  km  of  the  ship! 


7  Conclusions 

The  problem  of  association  and  fusion  of  radar  and  IRST  sensor  reports  is  not  straight¬ 
forward,  especially  because  IRST  provides  only  angular  measurements  while  radar  provides 
range  and  range-rate  measurements  as  well.  In  this  report,  simulation  results  have  shown 
that  a  centralised  EKF  tracker  is  a  solution  that  can  capitalise  on  the  higher  angular 
accuracy  of  the  IRST  sensor  to  provide  improved  track  accuracy  performance. 
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Appendix  A  Formulae  for  Converting 
Measurements  From  Spherical  to  Rectangular 

Coordinates 

The  formulae  are  due  to  Suchomski  [1999].  Denote  range,  azimuth  and  elevation 
measurements  by  r,  a  and  e  respectively.  Denote  measurement  error  standard  deviation 
by  oy,  a a  and  ae.  Then  the  converted  measurement  (x,y,z)  and  measurement  covariance 
R  are  given  by: 

x  —  r  cos  a  cos  e[exp(— cr2  —  cr2)  —  exp(—cr2/2  —  cr2/2)]  (Al) 

y  —  r  sin  a  cose[exp(— cr2  —  cr2)  —  exp(— cr2/2  —  cr2/2)]  (A2) 

z  =  r  sine[exp(— (jg)  —  exp(— crg/2)]  (A3) 

Rxx  Rxy  Rxz 

R  =  Rxy  -Ryy  Ryz  (A4) 

_Rxz  Ryz  Rzz_ 


where 

Rxx  — 

\r  (fixfixy  &x 

aXy)  +  Vr^PxPxy  ~  OLxo.xy )]  exp(-2cj2  -  2a2) 

(A5) 

II 

5ft 

0$ 

[r2(Pxy  -  a-xy  exp  a2)  +  a2 (1 2f3xy  -  axy  exp  a2)]  sin  a  cos  a  exp(-4cr^ 

-  2<r2) 
(A6) 

Rxz  = 

[r2(l-exp<7e2) 

+  cr2(2  —  expcr2)]  cos  a  sin  e  cose  exp  (—cr2  —  4  cr2) 

(A7) 

Ryy  — 

[r  (PyPxy  ~  Oty 

&xy)  A  Cr{C2f3y /3Xy  OtyOtxy)\  6Xp(  2(Ja  2(Je) 

(A8) 

II 

*4 

5ft 

[r2(l  -exp<72) 

+  cr2(2  —  expcr2)]  sin  a  sin  e  cose  exp  (—cr 2  —  4  cr2) 

(A9) 

Rzz  = 

[r2((3z  -az)  + 

(2 Pz  -  OLz)]  exp(— 2cTg) 

(A10) 

where 

&X 

=  sin2  a  sinh  cr2  +  cos2  a  cosh  a2 

(All) 

ay 

=  sin2  a  cosh  cr2  +  cos2  a  sinh  cr2 

(A12) 

Oiz 

=  sin2  e  cosh  a2  +  cos2  e  sinh  cr2 

(A13) 

&xy 

=  sin2  e  sinh  cr2  +  cos2  e  cosh  cr2 

(A14) 

Px 

=  sin2  a  sinh(2cr2)  +  cos2  a  cosh(2cr2) 

(A15) 

Py 

=  sin2  a  cosh(2 cr2)  +  cos2  a  sinh(2cr2) 

(A16) 

Pz 

=  sin2  e  cosh(2cr2)  +  cos2  e  sinh(2cr2) 

(A17) 

Pxy 

=  sin2  e  sinh(2cr2)  +  cos2  e  cosh(2cr2) 

(A18) 
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